First experimental determination of the charge radius of 74 Rb and its application in 

tests of the unitarity of the CKM matrix 
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Collinear-laser spectroscopy with bunched- beams technique was used for the study of neutron 
deficient Rb isotopes, out to 74 Rb (TV = Z — 37) at TRIUMF. The measured hyperfine coupling 
constants of 76 ' 78m Rb were in agreement with literature values. The nuclear spin of 75 Rb was 
confirmed to be I — 3/2, and its hyperfine coupling constants were measured for the first time. 
The mean-square charge radius of 74 Rb was determined for the first time. This result has improved 
the isospin symmetry breaking correction term used to calculate the Tt value, with implications for 
tests of the unitarity of the Cabibbo-Kobayashi-Maskawa matrix. 



The unitarity of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix is one of the cornerstones of the Stan- 
dard Model of particles and fields. Therefore, the value 
of all its matrix elements should be measured as accu- 
rately and as precisely as possible for the search of new 
physics phenomena, such as the existence of right-handed 
currents, extra Z bosons or another quark generation 
[l[. Indeed, over the years, the value of these matrix 
elements have been independently determined from dif- 
ferent experiments. Superallowed Fermi transitions of 
nuclear beta decays (0 + — > + ) provide the most com- 
petitive method to determine the up-down quark matrix 
element V uc i Q ■ Assuming the Conserved Vector Current 
(CVC) hypothesis, which states that the vector coupling 
constant of semileptonic weak interactions Gy remains 
unchanged in all superallowed beta transitions, this ma- 
trix element can be obtained from V u d = Gy /Gf, where 
Gf is the weak interaction constant for leptonic muon 
decay, and Gy is defined in terms of the corrected Ft 
value of the beta decay 0] , 
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The ft value is related to the measured half-lives, branch- 
ing ratios and the Q value of the beta decay. The terms 
S' R and 5ns are transition-dependent theoretical correc- 
tions, and 8c is the isospin symmetry breaking term 
The latter is purely nuclear structure dependent, the er- 
ror of which normally dominates the uncertainty on the 
Tt value. The term is transition-independent and 
K is a constant. Evaluations of the isospin symmetry 
breaking correction involve the 13 most precisely mea- 
sured superallowed beta emitters which span a wide Z 
range (6 < Z < 37) [4]. Within the existing theoretical 
approaches used to determine 5c [E[ the leading source 
of uncertainty stems from the incomplete knowledge of 



the radial extension of the charge distribution of the par- 
ent and daughter isotopes. Where experimental data are 
unavailable, the nuclear structure input is based on ex- 
trapolations of the existing models. For the lighter ele- 
ments, shell model calculations with a Saxon- Woods po- 
tential provide a more consistent framework from which 
this information can be reliably extracted. However, with 
higher Z , the orbitals of the fp shell are filled and shell 
model calculations become computationally difficult, due 
to the larger model space required. In addition, effective 
interactions become less reliable Q. There have been 
many theoretical efforts to fine-tune the shell model cal- 
culations in the fp shell, and the importance of the c/9/2 
orbital in this region has been acknowledged Q. Only 
by testing and refining these calculations against exper- 
imental data can a theoretical framework be developed 
for use in the calculation of 5c- In global evaluations 
of the Ft value d, the 74 Rb {N = Z = 37) isotope is 
the heaviest case included. It has the biggest 5c value 
among the existing cases (due to the Z dependence) , and 
the uncertainty around the 5c dominates the total un- 
certainty of the Ft value for 74 Rb at approximately 30% 
when all the fractional uncertainties are combined [l[ . 

The isotope shift of the hyperfine atomic spectrum is 
related to changes in the mean square charge radius [j| . 
The hyperfine structure allows one to extract the mag- 
netic dipole moment of the nucleus, fj,, if the nucleus has 
non-zero nuclear spin and if there exists a calibration 
isotope with a known moment. If the spin is I > 1/2 
the spectroscopic electric quadrupole moment, Q s can in 
principle be measured. Furthermore, the hyperfine struc- 
ture allows for a true measurement of the nuclear spin. 
High resolution laser spectroscopy is a powerful tool for 
probing the hyperfine atomic spectra, and the collinear 
method with fast beams has the sensitivity and resolu- 
tion required for radioisotopes Q. On the other hand, 
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low production (~ 10 4 ions/s) and background have be- 
come limiting factors for the sensitivity of this type of 
experiments. The use of cooler-buncher traps in con- 
junction with the collinear method is an elegant solution 
to this problem, as demonstrated at the IGISOL 0, 
and ISOLDE [10] facilities. The ISAC facility houses 
TRIUMF's Ion Trap for Atomic and Nuclear science (TI- 
TAN) 11]. This setup comprises a series of ion trapping 
devices used in particular for high precision mass mea- 
surements and decay spectroscopy studies. The first of 
these devices is a segmented linear Paul trap which is 
placed vertically in the ISAC beam line 12(. This de- 
vice can deliver cooled and bunched beams in forward, 
and uniquely, in reverse mode. The successful integration 
of the reverse-extraction of radioactive bunched ions for 
the collinear laser spectroscopy program has been demon- 
strated in the facility and is described in detail in [l]| ■ In 
this Letter, we present the first application of this tech- 
nique for the determination of the charge radius of 74 Rb. 
This result constitutes the first experimental benchmark 
of the theoretically determined 5c- 

Radioactive rubidium atoms were produced at the 
ISAC facility from spallation of a stack of natural Nb 
foils (density p ~ 8.5 g/cm 3 ) irradiated by 98 pA of pro- 
tons at 500 MeV from TRIUMF's cyclotron. The Rb 
atoms were surface-ionized, extracted at 28 keV, mass- 
separated and sent to the main experimental hall. The 
yield of 74 Rb measured at a tape station was 1.7 x 10 4 
ions/s. The ions were injected into TITAN's cooler- 
buncher and subsequently reverse-extracted to the laser 
spectroscopy setup, as described in reference (ljj . The 
ion beam was neutralized in a charge- exchange cell filled 
with Na vapor heated to 445 °C (neutralization efficiency 
of 50%). A frequency stabilized diode laser tuned to the 



of the ion bunch. The background suppression achieved 
was four orders of magnitude. 

The hypcrfine structure of both the ground state (I = 
0, Ti/2 — 17.66 m) and the isomeric state (E = 103 keV, 
7 = 4, T1/2 = 5.74 m) of 78 Rb are shown at the top of 
Figure [1] (a) . The single peak at the center of the figure 
corresponds to the ground state spin of 78 Rb (I — 0, and 
hence exhibits no hyperfine splitting). A \ 2 minimiza- 
tion routine was used for fitting the spectra using a Voigt 
profile. The same figure shows all the six hyperfine tran- 
sitions of the isomeric state. For the hyperfine structure 
of the isomer, the background and lineshape were fixed 
to the values obtained from the fit of the ground state 
transition. The plot also shows the hyperfine structure 
of 76 Rb together with the fit. As for 75 Rb, its measured 
hyperfine structure permitted an unambiguous ground- 
state spin assignment. Previously, the ground-state spin 
for this isotope had only been tentatively assigned as ei- 
ther I = 3/2 or I = 5/2 based on beta and gamma decay 



studies 15, [l6j]. A comparison of the \ 2 between the 
two assumptions confirms the 7 = 3/2 assignment. The 
bottom of Figure [TJ (a) shows the fluorescence spectrum 
of 74 Rb in a 400 MHz scanning region. We have excluded 
the resonance being in the region following the trend of 
the heavier isotopes. Since this isotope has ground-state 
spin I — 0, only one peak is observed. The signal-to- 
noise ratio of this peak is 2.2. The fitted FWHM is 
60(17) MHz, which agrees with the systematics of the 
measured widths of all the measured hyperfine transi- 
tions of 75,76,78m,78 Rb The regu i ts in Ta bi e [j] summa- 
rize the hyperfine constants and the isotope shifts of the 
Rb isotopes measured in this work. The hyperfine cou- 
pling constants obtained for 76,78m Rb, and the isotope 
shifts relative to 78 Rb confirm the previously measured 
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5s S'1/2 —> 5p -P3/2 transition, the U D 2 line" (780 nm) results reported in [14|. The isotope shifts 8v AA were 



was used, with 1.5 mW of light coupled to the beam line 
via an optical fiber, focused to approximately 1.5 mm di- 
ameter, and overlapped the atomic beam in a collinear 
geometry. For this purpose, two retractable apertures 
(1 mm and 3 mm diameter) separated by approximately 
12.5 cm were introduced in the detection region. The fre- 
quency locking of the diode laser was based on a polariza- 
tion stabilized He-Ne laser, using the principle described 
in reference [13|. The ion beam was decelerated by a 
voltage applied to the charge-exchange cell and Doppler- 
shiftcd into resonance. The laser induced fluorescence 
at 90° was recorded with a red-sensitive photomultiplier 
tube. An interference filter was used to help reduce the 
detected laser scatter. The cooler-buncher trap operated 
with a repetition rate of 10 Hz. Within this 100 ms cycle, 
80 ms was used to accumulate the ions into the trap and 
1 ms was used for cooling. The trap was left open for 
100 p,s for releasing the ions, and the extracted bunches 
had a temporal length of 4 ps at its full width at half 
maximum (FWHM). A 12 ps hardware gate was set on 
the detected photon signal in coincidence with the arrival 



TABLE I: Summary of the isotope shifts and hyperfine cou- 
pling constants (in MHz) of the Rb isotopes measured in the 
5s 2 Si/2 — > 5p 2 Ps/2 transition. For 75 Rb, the results of the 
hyperfine coupling constants obtained from the fits allowing 
for the two spin possibilities are shown. 
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78 - - - 0.83 

78m 4 +69.4(21) +1185.1(05) +29.3(01) +83.1(22) 0.76 

76 1 -24.3(12) -693(08) -17.15(02) +32(07) 1.08 

75 3/2 -41.1(17) +719.6(10) +17.8(01) +63(27) 1.18 

5/2 +12.1(12) +478.1(10) +11.9(02) +85.5(22) 1.87 

74 +99(10) - 0.93 

determined from the relative differences of the centroid 
frequencies between the isotopes (A) and the isotope of 
reference (A r ), 78 Rb. Before and during the data col- 
lection of 74,75,76 Rb, calibration measurements of 78 Rb 
were performed and the system has shown to be stable 
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FIG. 1: (Color on-line) (a) Optical spectra and fits for the neutron deficient Rb isotopes, measured on the D2 line, relative to 
78 Rb. The insert is an expanded plot of the 74 Rb dataset. (b) Isotope shifts of rubidium isotopes relative to A = 87. The 

,78m- ~~ 



black circles indicate the dataset presented in 
isotope shifts measured in this work. 



[141 ]. with the black asterisk indicating 5v 



-87 



Rb. The blue triangles are the 



to 1 MHz/day. Figure [U(b) shows the isotope shifts of 
rubidium isotopes relative to A — 87. It is interesting 
to point out the significant deviation of 74 Rb from the 
trend. 

The theory of optical isotope shifts separates the ob- 
served change in the transition frequency between two 
different isotopes with atomic masses A and A 1 into the 
mass shift and field shift terms, 



Sv 



AA 



5v AA 



= K 



MS ' 

A' 



A 



AA 1 



F5(rl h 



A A 



(2) 



The mass shift constant K is separable into the nor- 
mal and the specific mass shift. The first is a reduced 
mass effect and is exactly calculable, whereas the sec- 
ond is an effect of the change in the correlations be- 
tween pairs of electrons for the different isotopes. The 
quantity F is the field shift constant and depends on 
the atomic structure. In the case of the D 2 transition 
of rubidium, recent calculations performed by Dzuba- 
Johnson-Safronova produce F = —567.45 MHz/fm 2 with 



1% uncertainty [17|. Since the isotope shifts estab- 



lishes a relationship with changes in the mean square 
charge radius, the charge radius of a usually stable ref- 



erence needs to be used in order to extract an absolute 
value for the radius of the isotope of interest. The rel- 
ative rms charge radii of 85,87 Rb have been measured 
by muonic X-rays to be 85 Rb : (r 2 ^) 1 / 2 = 4.2035(3) 
fm and 87 Rb (reference isotope) : (r^) 1 / 2 = 4.1985 fm 
[3, such that 5(r 2 h ) 8r °- 87 = 0.042(2) fm 2 . With the in- 
formation that 6v S5 - S7 = -77.992(20) MHz at the D 2 
transition the mass shift constant on Equation [2] 

is K = -200(12) GHz/amu. The work of Thibault et 
al. gives 8v 7s ~ 87 = -478.4(15) MHz at the D 2 tran- 
sition Therefore, it is straightforward to tie this 
result together with our finding presented in Table U in 
order to obtain the isotope shift between 74 Rb and 87 Rb, 
viz. <^ 74 ~ 87 = -379(10) MHz. Inserting this value 
into Equation [H gives $(r 2 h } 74 ' 87 = -0.04(5) fm 2 , and 
therefore the rms charge radius of 74 Rb is determined 
to be (r 2 ^) 1 / 2 = 4.19(1) fm. The need for microscopic 
calculations to explain this result is motivated by the 
fact that admixtures of two V = + configurations with 
prolate/oblate deformations have been observed in the 
daughter nucleus 74 Kr. The strength of the E0 transition 
between these states was used as input to quantitatively 
understand how two-state admixtures are manifest in the 
charge radius [2(| . The observation of the isomeric decay 
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of 74 Kr by Chandler et al. has led the authors to sug- 
gest that the ground state wavefunction is composed of 
admixtures of states of opposite quadrupole deformation 
/3 2 = —0.35 and fi 2 — +0.38. The authors suggested that 
the mixing amplitude should lie between 30 — 50% [2l[ . 
On the other hand, there has been so far no experimen- 
tal evidence in support of a competing + state in 74 Rb 
22 , 23| . Shape coexistence and proton- neutron (pn) pair- 
ing correlations are expected to play a role in the JVftiZ 
region 2J], and the isospin dependence of the pairing 
interaction is likely to affect the excitation energy of the 
first excited + state. Since the charge radius is sensi- 
tive to static and dynamic deformation effects, the result 
for 74 Rb reported in this work can be used to assess the 
composition of the ground-state wavefunction. 

The main impact of the measurement presented in this 
work lies in the reevaluation of the isospin symmetry cor- 
rection term. In the analysis performed by Towner and 
Hardy 0] the authors split 5c into two contributions, 
5c = &C\ + Sc2- The first term 5ci reflects differences 
in isospin mixing between the parent and daughter iso- 
topes, whereas the second term 5c2 accounts for the ra- 
dial overlap correction between parent and daughter iso- 
topes. Although Sci is sensitive to the charge depen- 
dence of the effective interaction, its contribution is ap- 
proximately 10% of Sc2- The published value of Sc2 for 
74 Rb Q uses an extrapolation for the rms charge radius 
(( r 2ji/2 = 4.18(10) fm). This value was then used to 
fix the radius parameter in the Saxon- Woods central po- 
tential in order to calculate the radial-mismatch factor. 
For the computation of 5q2 the error budget comprises 
three pieces - (1) the component related to the mean 
square charge radius, 0.15; (2) the spread from Meth- 
ods II, III and IV, 0.20 (see ref. 0) and (3) the spread 
among the effective interactions considered, 0.12. Com- 
bining in quadrature and rounding the result gives the 
error of 0.30, viz. S C 2 = 1.50(30)%. With the mea- 
sured mean square charge radius of 74 Rb presented in 
this work, the first error reduces by an order of magni- 
tude, while the other two error components remain un- 
changed. This results in a reduction of the error on 5c2 
by ~ 20%, and hence reduces the uncertainty on the 
corresponding 5 c value. When one considers the error 
budget for the Tt value, this result puts the uncertainty 
on 5c to levels comparable to the uncertainty of the Q- 
value, which is the next leading term 25 -2l]]. Therefore, 
our result stimulates further experimental work to reduce 
the uncertainty of the latter. Indeed, this work has moti- 
vated recent precision Penning trap mass measurements 
on highly charged 74 Rb performed at TRIUMF. 

In summary, this work reports on the first measure- 
ment of extremely neutron deficient rubidium isotopes 
by optical methods. A novel experimental approach was 
developed and used to achieve unprecedented levels of 
sensitivity in the measurements of short-lived and low in- 



tensity radioisotopes produced at IS AC. The information 
obtained on the charge radius of 74 Rb was used to reduce 
the uncertainty in the evaluation of the isospin symmetry 
breaking correction term. This is the first time a charge 
radius measurement is being used for validating this cor- 
rection, in what constitutes one step towards improving 
the precision and accuracy of V u d , for tests of the unitar- 
ity of the CKM matrix. This experimental program will 
now be extended to measure the charge radius of 62 Ga, 
the second heaviest superallowed beta emitter. 
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